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We study static properties of attractively interacting Bose-Fermi mixtures of uniform atomic 
gases at zero temperature. Using Green's function formalism we calculate boson-fermion scattering 
amplitude and fermion self-energy in the medium to lowest order of the hole line expansion. We 
study ground state energy and pressure as functions of the scattering length for a few values of 
the boson-fermion mass ratio mb/mf and the number ratio Nb/Nf. We find that the attractive 
contribution to energy is greatly enhanced for small values of the mass ratio. We study the role 
of the Bose-Fermi pair correlations in the mixture by calculating the pole of the boson-fermion 
scattering amplitude in the medium. The pole shows a standard quasiparticle dispersion for a Bose- 
Fermi pair, for mb/mf > 1. For small values of the mass ratio, on the other hand, a Bose-Fermi 
pair with a finite center-of-mass momentum experiences a strong attraction, implying large medium 
effects. In addition, we also study the fermion dispersion relation. We find two dispersion branches 
with the possibility of the avoided crossings. This strongly depends on the number rario Nb/Nf. 



INTRODUCTION 



Recent developments in the field of cold atomic gases 
have proven that this system provides an ideal labora- 
tory for the studies of quantum many-body systems pj. 
This is due to the experimental facilities which allow to 
control various parameters characterizing the many-body 
system, e.g., external potentials including optical lattices, 
choice of atoms obeying Bose or Fermi statistics and their 
mixtures, variable particle densities, etc. The use of Fes- 
hbach resonances, in particular, makes the atomic gases 
an extremely flexible system as it provides a means to 
control atomic interactions^ l3(- One thus was, for in- 
stance, able to study the BEC-BCS crossover process in 
the two-component Fermi system, which has been under 
intense investigation for decades By changing the res- 
onance energies through the external magnetic field, one 
can in principle change the magnitude and the sign of 
the scattering length of the interacting particles, keep- 
ing track all the way from a resonating fermion pair to a 
bound composite particle, a bosonic molecule. 

The aim of the present paper is the study of pair corre- 
lations in a different system, a Bose-Fermi(BF) mixture. 
Degenerate mixtures of bosons and fermions have been 
created since several years, and studies of static and dy- 
namic properties have been performed^. Among those 
are the studies of attractively interacting BF systems, 
where one finds a sudden loss of fermions as the BF at- 
tractive interaction is effectively increased |(|. Detailed 
studies of the dynamics of this system are still missing, 
however. Recently the finding of Feshbach resonances 
and formation of the boson-fermion molecules have been 
reported [7|,|8j]. It is thus expected that by controling the 
BF interaction one may realize an analog of the process 
found in two-component Fermi systems. What should 



be expected if one replaced fermion pairs in the BEC- 
BCS crossover process by BF pairs? Such studies have 
indeed been performed theoretically [i| (see also By 
adopting a Cooper type two-particle problem on top of 
the boson-fermion degenerate system, it was shown that 
a stable correlated BF pair is created even before the 
threshold for the BF bound state. In contrast to the BCS 
case, however, the system allows only one correlated BF 
pair with a given center-of-mass (CM) momentum be- 
cause of the fermionic nature of the composite particle. 
It is then suggested that by increasing the BF attractive 
interaction, one may create BF pairs with different CM 
momentum stepwise, until finally a new Fermi sea of the 
BF pairs is completed. 

In Ref. [9( a separable BF interaction has been adopted 
to elucidate the mechanism of the creation of BF pairs. 
In the present paper we adopt a standard pseudopoten- 
tial for the interaction, and calculate energy and pressure 
of the system for various values of input parameters. We 
use Green's function formalism for this system and calcu- 
late perturbatively relevant diagrams to lowest order of 
the hole-line expansion. Such formalism has been devel- 
oped in [TTj] together with the calculation of the energies 
including Bosc-Bose (BB) interaction. Our formulation 
is similar to [Til ], but we use the renormalization proce- 
dure of [TH in relating the pseudopotential strength to 
the iS-wave scattering length. This allows us to formally 
take the limit |a| — > oo, the unitarity limit [T3j . which 
is necessary when one considers a (nearly) bound state 
of a pair of atoms. We then calculate the poles of the 
BF pair scattering amplitude in the BF medium, which 
may be compared with the results of @ . Studies of the 
behavior of the poles as a function of input parameters 
give us suggestions on the role of the BF pair correlations 
in the static properties of the system. 



2 



The content of the paper is as follows: In the next 
section we present our model based on the Hamiltonian 
without Bose-Bose interaction. We calculate the BF scat- 
tering amplitude in the BF mixture in ladder approxima- 
tion, and give formulas for physical quantities in terms of 
the amplitude. In section 3 we show numerical results for 
the ground state energy and pressure for various choices 
of the boson/fermion masses and the values of the Bose- 
Fermi interaction. We then study Bose-Fermi pair cor- 
relation in Section 4 by focusing on the pole structure of 
the boson-fermion scattering amplitude in the mixture. 
We also calculate the pole of the single fermion Green's 
function and study the role of the Bose-Fermi pair and 
its dispersion. We summarize our results in section 5 to- 
gether with a comment on the effects of the Bose-Bose 
interaction. Detailed expressions for the scattering am- 
plitude are given in the appendix. 



II. FORMULATION 



omit the fluctuation of the boson number in the conden- 
sate. The boson number operator iVf, writes 



^6 = ^0+ / dW(*Mx). 



and the Hamiltonian takes the form 
H = H + Htf , 



where 
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We consider a uniform system of a polarized Bose- 
Fermi mixture of atomic gases with attractive boson- 
fermion interaction. The model Hamiltonian of the sys- 
tem is given by 



H = T b + Tf + H b f , 



T, 



S 



d a x^ + (x 



V 2 
2nib 
V 2 
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MfcJ 0(x), 
V>(x), 



Hbf = 9bf / d 3 x0t(x)^(x)VXx)0(x), 



(1) 



where tp and </> are the boson and fermion field operators, 
respectively, Tf, and Tf denote bosonic and fermionic ki- 
netic energies, while Hbf denotes boson-fermion interac- 
tion with strength gbf (< 0) of the boson-fermion pseu- 
dopotential. Effects of the boson-boson interaction will 
be mentioned later, while the fermion-fermion interaction 
is omitted throughout as we consider one-component (po- 
larized) fermions. We will adopt the Bogoliubov approx- 
imation in treating the Bose-Einstein condensate (BEC), 
and therefore include in Tf, the bosonic chemical potential 

A*6- 



H b f = n Q g b f j d 3 X7/; t (x)i/'(x) 

+ Vno9bf J d 3 x^ t (x)V'(x)((^ t (x) + <p(x)) 

+ g bf fdV(x)^(xMx)^(x). (6) 



Physical quantities can be expressed in terms of Green's 
functions. We define the boson and fermion Green's func- 
tions by 



iG* {x — y) = 



iG\x - y) 
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(7) 
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where ipn (x) , ips (x) are the field operators in the Heisen- 
berg picture, and |\&o) represents the interacting ground 
state. 

The energy of the system can be expressed in terms of 
the Green's functions. The fermion and boson kinetic en- 
ergies are calculated according to the standard procedure 
[IJI as: 



A. Green's function formalism in the Bose-Fermi 
mixture 

To treat condensed bosons, we adopt the conventional 
Bogoliubov method by separating the zero momentum 
mode from the remainder : 



(2) 



770 + ¥>(x) 



together with its conjugate. Uq = N /V is the number 
density of bosons with momentum k = 0. As usual we 
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IQOV^ ( 9 ) 
'' r'bf „\„iqan 



where G(q)'s are the Fourier transform of the Green's 
functions, 77 is a positive infinitesimal, and we set e b ^ = 
q 2 /277if,j. The different signs in the two expressions 
come from the ordering of the field operators. 

To calculate the interaction energy, we first consider 
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the Heisenberg equation of motion for the fermion field: 

d 

i^feW =Wh(x),H] 
V 2 

= (- o )^h{x) + n g b ft/jH(x) 

2m/ 

+ y/no9bf1pH(x)((p]j ( x ) + ^ffW) 

+ 9bf1pH(x)(PH(x)(p H (x). (11) 

Multiplying by ip H (x') and integrating over x, we obtain 



(H 



bfi 



• fd 3 x lim (i^- + — ^ G f (tx,t'x') 
J x'^x V 9t 2m f J v ; 



-iV 



d 4 g 
(2^ 



(12) 



We now introduce fermion and boson self-energies £^(<?) 
and S 6 ((?) through 



and 



G f {q) = p , 



G"(g) = 



g -< + Mfc-£ b (g)' 



(13) 



(14) 



In the integrand of Eq. (|12[) one may use the relation from 
Eq.CSJ 

( 9o - e £)G%) = l + E/ (<?), (15) 



and finds 



= -iV 



d 4 q 



V f (q)G f (q)e iqoT >, (16) 



the first term in the r.h.s. in Eq. (|15p giving a null contri- 
bution to the integral. The total energy E of the system 
is finally obtained as 

E =(7» + (T b ) + (H bf ) 
d*q 



iV 



iV 



(2tt) 4 
d 4 q 
(27T) 4 ' 



{e{ + &(q))Gt(q)e i «"' 
b cl G b (q)e lqorl . 



(17) 



The thermodynamic potential at zero temperature is 
given by 



Sl(N f ,N ,n b ) = {H)=E- f i b (N b ), (18) 



where 



(N b ) =N Q + l vj J^L G \q)e^. (19) 

The system is characterized by the boson and fermion 
particle numbers, N b and Nf. The number of bosons 
satisfies the thermodynamic relation 



d(J, b 



-N b 



(20) 



The parameter Nq should be chosen to minimize the ther- 
modynamic potential 



on 

dN~ 



(21) 



which leads to an explicit expression for fi b as shown 
below. 

We also will calculate the pressure to discuss the sta- 
bility of the system: As usual, it is obtained from the 
thermodynamic relation 



P 



dE 

w 



(22) 



B. Self-energy in the ladder approximation 

To obtain the total energy of the system we calcu- 
late the fermion self-energy T, F in ladder approximation. 
Here the self-energy is expressed in terms of the two- 
particle scattering amplitude, T(q, q', P), in the medium 
of the Bose-Fermi mixture as shown in figHJ The inter- 
action energy is accordingly calculated up to the lowest 
two-particle correlation diagram, figlU in the spirit of the 
hole-line expansion [l4l |. 
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FIG. 1: Fermion self-energy in ladder approximation. The 
double solid lines represent fermion propagation, while single 
solid line represents noncondensed free boson propagation. 
The arrows denote condensed bosons and are associated with 
the factor y/no. The zigzag lines represent the boson-fermion 
interaction gtf. 
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FIG. 2: Ladder diagram contribution to the interaction en- 
ergy. The downward double solid line indicates a hole propa- 
gation. Otherwise as in fig[T] 

The scattering amplitude T in the present model obeys 
the integral equation [§, [ll], [lj, [ll| , 
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r(q,q',-P) = gbf + igbf 



(2ir)^° 



mf + m b 



P + k)G° 



r 



m b 



i f + m b 



P-k T(k,q',P), 



(23) 



where P denotes a four-momentum of the center-of-mass 
motion of the interacting particles, while q and q' are the 
relative three-momentum in the final and intial states. 
The boson and fermion free Green's functions in medium 
are given by 



G f Q (p) 



e(\p\-k F ) , 6(k F -\ P \ 



Po 



J 



+ if] Po 



ii] 



(24) 



G b (p) 



Po 



]i b + it] ' 



(25) 



where the Fermi momentum kp is fixed by the fermion 
density Nf/V. After the integration over k , Eq. 
becomes 



r(q,q',P) =g bf +g bf 



d 3 k e(\P f + k 



(2^) 3 p -4 



p/+k p b -k 
i 



— r(k, q ',p) 



(26) 



with P/ = m,f/(mf + m b )V and P b — m b /(mf + m b )P. 
We dropped the hole propagation part in accordance 
with the present approximation. With respect to the 
T- matrix equation in [§], we notice that there the phase 

space factor in (f2"6"j) is replaced by 8 ( P/ +k — kp^j ^ 

9 ^ Pf + k — kFj + No- This is natural, because in [§| 

the shift operation © for the bosons has not been per- 
formed and therefore the free boson occupancy ap- 
pears additionally. The two formulations are, however, 
essentially equivalent. From the structure of Eq. ([26|) . one 
easily finds that T depends only on the variable P, and 
we hereafter write simply r(P). One also finds that the 
integral in Eq. (|26p requires a momentum cutoff, which 
originates from the use of the zero-range interaction. We 
can remedy this shortcoming by employing the observ- 
able S-wave scattering length a, instead of the pseudopo- 
tential coupling constant g b / . We perform this renormal- 
ization following the procedure adopted in [12| (see also, 
[ill), slightly different from the one in The S-wave 
scattering length is related to the two-particle scattering 
amplitude Tq in vacuum by the relation 



To(q 



P = 0) 



2na 

3 

V 



where v is the reduced mass, and Tq obeys the equation 
similar to Eq. ([23]) with Go replaced with the free Green's 
function in vacuum. By solving the equation for To one 
obtains 



2na 



9bf 



9bf 



d 3 k 



(28) 



(27T) 



3 J 



where the integral in the denominator involves again the 
implicit momentum cutoff. Now one may combine the 
above expression with Eq. (|26p . and eliminate g b f in favor 
of the scattering length a, and finally obtains 



r(P) = 



2na 



l-^I(Po,|P|) 



(29) 



(27) with 



I{P^\P\) 



(2^ 



P/+k 



+ 



4 + 4 



(30) 



Since the integral in the denominator is convergent at large |k|, we can let the momentum cutoff go to infinity. 
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The expression (|29[) involves all orders in the scatter- 
ing length and allows us to formally take the unitarity 
limit \a\ — > oo in the following section. This limit has 
been studied for two-component Fermi systems in rela- 
tion with the BEC to BCS crossover phenomenon. If a 
similar phenomenon is expected or not for Bose-Fcrmi 
pairs will be studied in the next section. 

Using above vertex function, we can calculate the 
proper self-energies for the fermion and the boson as 

&(p)=n r(p) (31) 

Z b (p) = -ij 0^G f o (p')T(p + p'). (32) 

Expression (f2"6"| implies that T(p + p') is analytic in the 
upper half p' plane, and Eq. (f3"2")) reduces to 

z b ( P ) = J ^e{k F - \p'\)r(p+ P '). (33) 

with p' = Cp,. This shows that £ b (p), and hence, also 
G b (p) is analytic in the upper half po plane. One then 
finds from Eq.(20|) that 

N b = N . (34) 

III. RESULTS FOR TOTAL ENERGY AND 
PRESSURE 

We calculate the energy of the system in the leading or- 
der of the hole-line expansion, that is we replace Green's 
functions in Eq. (fT7|) with the free one Gq in Eq. ([23|) . 
and obtain 

E~-ivjj0^ + E'(«)) Gg(g)e 4 «"» 

- E N f + N J ^0{k F - |p|)r(e£, p), (35) 

where, Eq — 3/5Ep and Ep is the Fermi energy. Within 
the same approximations, the thermodynamic potential 
at zero temperature is given by 

n = E N f +N / -^0(fcF-|p|)r(e£,p)-M6iVo- (36) 
J (27r) 

Thus, the equilibrium condition (12~1]) for fl leads to the 
integral equation for /if,, 

M6 = / ^^-|p|)r(e^,p), (37) 

where T depends also on fi b . The total energy of the 
system is then finally given by 

E = E N f + n b N b . (38) 



Details of the calculation and the analytic expression for 
r are given in the appendix. 

We may rewrite Eq. ((37|) in a scaled form as 

h = 2 (l + J dpp 2 f(p, fi b , C), (39) 

where we introduced tilde (dimensionless) quantities 
through r(p, /if,) = 2-n: 2 /vkp T(p, fib), a — kpa, p — 
|p|/fcp, fib = [ibjEp. The expression shows that the 
scaled chemical potential jib depends only on the mass 
ratio C = mb /m / and the dimensionless scattering length 
a = kpa. We solved Eq. (|39|) for fib numerically as a func- 
tion of the boson-fermion mass ratio C for different values 
of the interaction strength represented by a. In terms of 
the scaled quantities, the ground state energy per particle 
is expressed from Eq. ([38|) as 

A = | £f( i + /3)) (40) 
N f 5 

where the dimensionless parameter /3 is given by 

We first show the results for energy and pressure in the 
unitarity limit, \a\ — > oo. In this limit, assuming S-wave 
scattering and neglecting effective range, we are left with 
only one length scale, kZ 1 , or in terms of the den- 

sity nj — Nf/V [H, [TtJ for a given mass ratio C,. Note 
that the chemical potential fib has no kp dependence in 
the unitarity limit, and the parameter (3 depends only on 
the number ratio Nb/Nf. Thus the ground state energy 
per particle, Eq. (|40)l . is proportional to Ep, and the de- 
pendence on the parameters are all absorbed in a simple 
multiplicative factor (1 + (3). 
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m b /m f 

FIG. 3: Scaled chemical potential fib = fib/EF as a function 
of the boson-fermion mass ratio £ = mb/mf in the unitarity 
limit. 

We show in fig[3] the chemical potential Hb and in 
fig0] the beta parameter as functions of the mass ratio 
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0.5 1 1.5 



FIG. 4: f3 parameter as a function of the boson-fermion mass 
ratio in the unitarity limit for Nb/Nf = 2.0, 1.0, 0.5 from top 
to bottom. 



C = mb/mf, both in the unitarity limit. Note that the 
results are independent of the magnitude of the individ- 
ual mass parameters as we show dimensionless quantities 
scaled with Ep . Figure [3] shows that the boson chemical 
potential is always negative. This fact reflects the at- 
tractive boson-fermion interaction in the unitarity limit, 
in accordance with Eq. (|28[) which implies negative gbf. 
The behavior of in figH] simply follows the one of the 
chemical potential. The results suggest that the attrac- 
tive interaction becomes more effective for small values 
of the mass ratio £, and the effect is greatly enhanced 
as the particle number ratio Nb/Nf becomes larger, that 
is as the number of bosons increases with respect to the 
fermions. The dependence on £ may partly be under- 
stood by noting that the relative phase space available 
for the intermediate states in the two-body scattering in 
the mixture will be larger for small £, i.e., for a relatively 



larger m r , because of the lower Fermi energy and higher 
level density. 

We next consider the pressure to study the stability of 
the system. Since the total energy takes a universal form 
and the (3 parameter has no volume dependence in the 
unitarity limit, the pressure is simply given by 



P 



dE 
W 



2_Nf_ 
5 V 



E F (l + 0) (unitarity limit). (42) 



For (3 < — 1 the pressure becomes negative, and the sys- 
tem collapses. This happens especially for larger values 
of Nb/Nf, where the pressure becomes always negative 
irrespective of the mass ratio £. In actual experiments, 
e.g., for the 40 K- 87 Rb mixture which has £ ~ 2.3, the 
number ratio is typically 0(1) ~ 0(1O 3 ) and the sys- 
tem would collapse in the unitarity limit. This is not in 
contradiction to recent experimental results Q . 

Another feature in the unitarity limit seen from fig[3] 
and figlU is that the boson chemical potential and the 
ground state energy saturate when the mass ratio be- 
comes large. It is natural that the bosonic degree of 
freedom gets frozen and a universal fermionic description 
appears in this case, since the fermionic effects on the 
bosons would become negligible, and the bosons would 
act as a static external field for fermions. This behavior 
at large ( also has been discussed in [ll[ . 

We now study the case with an arbitrary value of the 
scattering length. First, we show the chemical potential 
fib as a function of (kpa)^ 1 in figO where the mass ratio 
C and Nb/Nf are set to 1. Then, we show the energy 
and pressure as a function of (kpa)^ 1 in figs O [7] and [51 
for the mass ratio £ = 0.8, 1.0, 1.2 with different values 
of Nb/Nf = 0.5, 1.0, 2.0. For a^0, the pressure is given 
by 



2^ 
5 V 



-E f 



2 1 



Nb 

N fJi; 



dpp 2 (f(p,fM b ,() 



^r 2 (p,/i h ,c) 

a 



(43) 



r 



where the term dependent on T 2 reflects that the chemi- 
cal potential, and hence the (3 parameter, depends on V . 
From the results on energies we see that the system be- 
comes more attractive at smaller values of the mass ratio 
£ as in the unitarity limit, although the effect is not large 
in this parameter range. We find a strong increase of the 
attraction as the parameter (kpa) -1 passes through zero, 
the unitarity limit, from negative to positive. This is in 
accord with a naive picture where the positive values of 
the scattering length imply a newly formed bound state. 
One should however note that even in the present approx- 
imation the effects of the medium modify the two-body 
scattering amplitude, and a simple picture of indepen- 



dent bound pairs does not hold in general. Turning now 
to the pressure, a comparison of Eq. (|43p with the cor- 
responding expression (|42p in the unitarity limit shows 
that the strong attraction for a > comes from the large 
negative values of T as well as from the coherence of the 
two terms in the integrand. 



IV. BOSE-FERMI PAIR CORRELATION 

Results of the previous section indicate that the strong 
attraction in the mixture will show up especially for pos- 
itive a, which may eventually lead to a collapse of the 
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FIG. 5: Scaled chemical potential fit = fib/Ep as a function 
of (/cp-a) -1 for uib/mf = Nb/Nf = 1. 



system. We now consider another scenario for the at- 
tractively interacting mixture, the possibility of a Bosc- 
Fermi pair formation [9j, [lCj. For this purpose we study 
in this section the behavior of the pole of the Bose-Fermi 
scattering amplitude T(P) in the mixture. 

In unpolarized (or two-component) Fermi systems, an 
infinitesimal attraction around the Fermi surface leads 
to formation of Cooper pairs with a center-of-mass (CM) 
momentum P = 0, causing a transition to the BCS state. 
In a Bose-Fermi mixture, on the other hand, the differ- 
ence in the momentum distribution of the two particles 
and the fermionic character of the Bose-Fermi pair, in 
particular, predict quite a different scenario for the for- 
mation of the pairs in the mixture. It requires consid- 
eration of the balance of the kinetic energies of different 
kinds of particles and the magnitude of the attractive 
interaction. 



A. Preliminary considerations 

Let us give a brief picture on the formation of Bose- 
Fermi pairs in the mixture. Following the idea of 0], we 
may take a Bose-Fermi mixture with Nb — Nf, focusing 
only one pair of a boson and a fermion with CM momen- 
tum P, putting other particles as a free background. The 
Hamiltonian in this system may be written as 



H 



Pi 



n 2 



2mt, 2m f 



Vbf + Hb g , (pb + Pf 



(44) 



where we explicitly write the kinetic energy and the inter- 
action for the two particles, while the background Hamil- 
tonian iJbg acts only to impose a Pauli-principle con- 
straint. We neglect here the effect of the boson chemical 
potential for simplicity. If the effect of Vbf on the pair 
were negligible, the energy of the Bose-Fermi pair with 
CM momentum P would be simply 



£free(|P|) = 



2m i 



(45) 



which has been called free branch in Q , since the boson 
will remain at Pb = in the condensate. When the effect 
of the interaction becomes important, one may rewrite 
Ea.pg]) as 



H 



2(m& + nif) 



Hrel + Hbg, 



(46) 



where the interaction Vbf is contained in the Hamilto- 
nian H re \ of the relative motion. By replacing the latter 
with its eigenvalue E re \ (the effect of the medium may be 
included here), one obtains a different dispersion curve 
for the Bose-Fermi pair, the collective branch, which is 
given by 



econ(|P|) 



2(m b + rrif) 



E, 



rel • 



(47) 



The calculation for the model with a separable interac- 
tion in Q shows that the two dispersion curves, Eqs. (l4"5l) 
and (f47|) , coexist except for the region of CM momentum 
P c given by a solution of the equation 



ef ree (|P c |) = e C oll(|Pc 



(48) 



There is a mixture of the two branches around |P| = 
P c , and the dispersion curve of the pair deviates from 
£f r eo and e co ii- For a sufficiently attractive interaction, 
the solution of Eq. (|4"8"|) satisfies the condition |P C | < kp, 
which occurs for 



■nib + mf 



-E F . 



(49) 



One may thus expect that the system will lower the en- 
ergy by converting the free Bose-Fermi particles in the 
range |P C | < |P| < kp into collective pairs. (Note that 
we neglect the interaction of the pairs in this simple ar- 
gument.) If Vbf is so strong as to allow for a bound state 
of the pair, i.e., E rc \ < 0, all the free bosons and fermions 
would be replaced with the bound B-F pairs, and the new 
Fermi sea of the pairs will be formed. 

This picture may be modified even in this simple 
model, however, for a system with many Bose-Fermi 
pairs. As the pairs can occupy low-momentum states 
having lower energies without changing the total momen- 
tum of the system, the formation of the collective pairs 
may start if the condition 

E lcl = e co ii(|P| = 0) < e hcc (k F ) = Ep (50) 

is satisfied, even before the condition (|4"9")) . Similarly, the 
formation of the Fermi sea of the pair will be completed 
only when 



ecoll(fcf) <£frcc(|P| =0)=0 



(51) 



is satisfied. 

We note that the collective branch (|47[) may suffer a 
Laudau type damping for sufficiently large E TC \. This is 
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FIG. 6: Scaled energy E/EpNf (left panel) and pressure PV / EfN f (right panel) as a function of (kFa) 1 for mi/rrif = 0.8. 
Three curves corresponds to Ni/Nf = 0.5, 1.0, 2.0 from top to bottom. The vertical line indicates the unitarity limit a — oo. 
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FIG. 7: Same as figlU but for mj/rn/ = 1.0. 



because the free Bose-Fermi mixture with Fermi momen- 
tum kp has a continuum of the boson-fcrmion excitation 
with momentum P = p& + p/ which starts at the energy 
1 

e th (\P\)= m ~ kF)2 +E F . (52) 
Zrrib 

This implies that the collective branch remains un- 
damped only when 

eoou(|P|) < e th (|P|), (53) 

a condition which is approximately realized for the pole of 
the boson-fermion scattering amplitude as shown below. 

B. Pole behavior of the two- particle scattering 
amplitude 

Now we study the behavior of the pole of r(P). A first 
study concerns the pole condition. That is 

^=/(Po,|P|). (54) 
lira 

For |P| = 0, we show in ng[5]the right hand side of (|54")) as 
a function of Pq for m^/mf = 1. We see the development 



of a logarithmic divergency as Pq aproaches k F /2v — [Ib- 
This stems from the fact that above dispersion integral 
has exactly the same structure as the one encountered in 
the problem of Cooper for a fermion pair in a Fermi-sea. 
We want, however, to point out that the pole corresponds 
to a composite fermion what has important consequences 
for the physics. Nevertheless the fact is there that a 
stable collective B-F pair developes for any infinitesimal 
attraction, i.e, even in the limit a — > — quite in analogy 
to the original Cooper pole 

Let us now discuss the collective pole contained in 
r(P). Be P C the pole of T(P) with CM momentum P. 
Pq(|P|) represents the total energy of a boson-fermion 
pair corresponding to the collective branch aside from the 
chemical potential. One may then define a |P|-dependent 
binding energy (including total kinetic energy) measured 
from the last filled free Bose-Fermi pair by 

A pQir (|P|) = e fre e(fc F ) - M6 - ^o c (|P|), (55) 

see figHUl Positive value of A pa j r (|P|) would signal a 
formation of the Bose-Fermi pair. 

We solved Pq (|P|) numerically. Here we include chem- 
ical potential fib-, and show A pa i r (|P|) in units of Ep 
for C = 0.8, 1.0, 1.2 in figsdU [TU and [H] The number 
ratio is fixed at Nb/Nf = 1. Left panel of each figure 
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FIG. 8: Same as figlU but for mt/rrif — 1.2. 




-2 -1.5 -1 -0.5 0.5 1 1.5 2 2.5 



Po 

FIG. 9: The behavior of the right hand side of eq.([54]). The 
horizontal line is the left hand side. Here, rrib/mf — 1. The 
vertical line shows the the position of the logarithmic diver- 
gence. 




FIG. 10: Rough sketch of the formation of the Bose-Fermi 
pair with CM momentum P. The dots represent a boson 
(left) and a fermion( right). The shaded box indicates a fermi 
sphere. 

shows results for several values of 1/kpa. One finds that 
each line extends up to a maximum value of |P|, the en- 
ergy of which corresponds to e t h(|P|) in Eq. ([5^|) . i.e., the 
point where the pole hits the continuum and obtains a fi- 
nite lifetime due to the non- vanishing imaginary part. In 
the right panels of these figures, we show the dispersion 
curves at the threshold values of 1/kpa, where A pa i,.(|P|) 



turns from negative to positive for the first time. 

We see from these figures that the effective binding 
increases as the interaction becomes more attractive, and 
eventually leads to a formation of the Bose-Fermi pair. 
The threshold value of is lower for a smaller value 

of £ = rrib/mf in agreement with the result for the total 
energy. We find also that A pa i r is always negative in the 
unitarity limit, suggesting that the Bose-Fermi pair may 
not be formed in this limit. 

A peculiar feature seen from the figures is the non- 
monotonic dependence of A pa i r against CM momentum 
|P|, which is in contrast to the P 2 dependence of Eq. (l47|) 
in the simple picture. This is particularly apparent for 
small values of £ = mt /m f . Missing in the simple model 
is the P-dependence of the pair binding energy E rc \ which 
should be present due, e.g., to the phase space available 
for the two-body scattering in medium. The medium 
effect becomes stronger for larger (positive) values of 
(fcpci) -1 , especially for small £ as seen from the figure. It 
may be mentioned in addition, that the results in the next 
subsection suggest a strong interplay of the collective and 
free branches for the Fermion dispersion relation. 



C. Fermion dispersion and level crossing 

Let us now investigate the pole structure of the single 
particle Green's function G$ of (fT5|) . The pole condition 
reads, 

£p = e£ + £'(£ p ,p). (56) 

In order to be consistent with our publication in 0], 
we here neglect the chemical potential Hb in Y(P) of 
(2HJ). This means that we replace in Y(P) the Boson 
and Fermion propagators by totally free ones with the 
kinetic energies of boson and fermion and ji^ = which 
is the uncorrelated value. 

With the definition of in (fT3|) with /!& = and the 
expression Y of P9")) , we easily can investigate the poles 
of G* for various system parameters. We first consider 
cases, like in Q, where Nb <C Nr. 
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FIG. 11: Apair(P) as a function of center of mass momentum |P| for various (fc^a) -1 with mj/m/ = 0.8 and Nb/Nf — 1. The 
shaded area shows the contiuum of the excitation which starts at e t h, Eg . (I52[l . The right panel shows the enlarged figure at 
the threshold value of (fcpa) — 1 (see text). Note the scale difference in the vertical axis between left and right panels. 




FIG. 12: Same as FiglTTlwith m b /m f = 1.0. 



We expect two branches: one which corresponds in 
week coupling to E p ~ e£ and one which corresponds to 
the collective branch, i.e. the pole contained in T(p). 



In figHH we show on the left panel the case of Nb/Nf — 
1/1000 and on the right panel the case of N b /N f = 1/100. 
From top to botom, we have (fcpa) -1 = —0.2,0.2,0.5. 
Here mass ratio mb/mf = 1. The dotted area is the 
region where the imaginarypart of is different from 
zero. As in our ealier work Q, we see an avoided crossing 
of the two branches at some finite value of p/kp. We 
also see that the interaction between the two branches 
becomes stronger as Nb/Nf increases. 



In figHU we show the case (£;/a) _1 = 0.5 and Nb/Nf = 
1/10 for left panel and Nb/Nf = 1 for right panel. We 
see that no crossing feature is visible any longer and the 
two branches probably become completely hybridised. In 
a future publication we intend to investigate in a system- 
atic way the nature of the two branches as a function of 
the system parameters. 



V. SUMMARY, DISCUSSION AND 
CONCLUSION 

We studied in this paper the static properties of the 
Bose-Fermi mixture in the lowest order of the hole-line 
expansion, i.e. in the T-matrix approximation. The in- 
teraction parameter is expressed in terms of the scatter- 
ing length up to infinite order using the renormalization 
procedure of ref.[T3|, so as to allow for the calculation 
around the unitarity limit. The T-matrix approach is a 
common approximation often used in the past in Fermi 
and Bose systems [lij]. It has, however, not been investi- 
gated very much in Bose-Fermi mixtures and, therefore, 
one has not much experience about its quality in that sit- 
uation. Recently appeared, however, a study for a one di- 
mensional system |18| , from where exact Quantum Monte 
Carlo results are available for comparison. In ref [l|| it 
is shown that the T-matrix approach yields quite reason- 
able results also in the Bose-Fermi case, even in a one 
dimensional case which is probably the worst situation 
possible. 

With this background in mind, we first studied energy 
and pressure as functions of the inverse scattering length 
for several choices of the mass ratio £ = mb/mf and 
the number ratio Nb/Nf. The energy of the system be- 
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FIG. 13: Same as FiglTTIwith m b /m f = 1.2. 



comes strongly attractive as the inverse scattering length 
changes sign from negative to positive, i.e. around the 
unitary limit. As one increases the number of bosons 
with respect to the fermions, arrives a point where the 
pressure becomes negative, i.e. the system becomes un- 
stable (collapse). The effect is stronger for small values 
of TJih Imr. This is not in contradiction with experiments 

1- 

Next we studied the possibility of stable BF-pairs, as 
in Q. We, indeed, also found in the present model that 
even for infinitesimal BF-attraction a stable BF-mode ap- 
pears, reminiscent of the Cooper pole in a two component 
Fermi gas, since in both cases its origin stems from the 
presence of a sharp Fermi surface. However, in the BF 
case the BF-pair is a composite fermion, whereas in the 
original Cooper problem, one has a composite fermion 
pair, i.e. a boson-like cluster. In the latter case, many 
pairs can well be treated by the usual BCS formalism. 
On the other hand , the case of many BF-pairs still has 
to be worked out. This shall be done in future work. But 
we here studied the BF pair formation by observing its 
binding energy measured from the last filled free BF pair. 
For some finite values of the attractive interaction, there 
occurs the formation of stable BF pairs. For ( > 1 the 
pair shows a standard dispersion of a quasiparticlc, while 
for £ = 0.8 the pair with finite center-of-mass momenta 
feel stronger attraction. This effect is not clearly seen 
in the energy or pressure, where the singular effect may 
have been averaged out. 

Up to this point we assumed an ideal case where 
the Bosc- Fermi interaction is dominant, while the Bose- 
Bose interaction was neglected. To compare with a real 
atomic gas system, the Bose-Bose interaction cannot be 
discarded even when the Bose-Fermi interaction is en- 
hanced, e.g., through Feshbach resonances. We checked 



the effect of the Bose-Bose interaction on the energy 

1/3 

of the system up to the first order in N b ' abb- We 
took the BoseBose scattering length a^b in the range 
—0.3 < kF<ibb _ 0.3, and the number ratio Nb/Nf =2.0, 
1.0 and 0.5, and repeated the calculation of the (3 parame- 
ter and the pressure 20] . The calculation shows that the 
effect of the Bose-Bose interaction is small within the 
adopted parameter values: The (3 value at the unitarity 
limit, for instance, deviates less than 10% for nib = m f 
and Nb — Nf, and does not change conclusions obtained 
at abb = 0. 

In summary, we suggest that the energy gain in the 
Bose-Fermi mixture at positive values of abf is related to 
the formation of the resonant Bose-Fermi pairs, and that 
the center-of-mass momenta of the pairs are dependent 
on the ratio nib /m / due to the statistics of the two kinds 
of the particles. 



Acknowledgments 

We thank X. Barillier-Pertuisel, K. Suziki, T. 
Nishimura, T. Maruyama and H. Yabu for useful dis- 
cussions. We also thank J. Dukelsky and S. Pittel for 
their general interest and contributions to the subject of 
BF correlations. 



APPENDIX: EXPRESSION FOR THE 
SCATTERING AMPLITUDE 

Let us calculate the I(Pq, |P|) of the scattering 
amplitude T(P) 



12 



LJJ 
LU 



2 
1.5 
1 

0.5 



- 






(k F a)" 1 = -0.2 








UJ 

" Q. 

UJ 



2 
1.5 
1 

0.5 







- 


(k F a)" 1 =-0.2 


I I I 






0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
p/k F 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

p/k F 



LJJ 
LU 










(k F a)" 1 =0.2 






0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
p/k F 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

p/k F 



LJJ 
LU 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
p/k F 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

p/k F 



FIG. 14: Fermion dispersion curves. Left row corresponds Nt/Nf = 0.001 and right row corresponds N^/Nf — 0.01. From top 
line to bottom, (fcFffl) -1 are —0.2,0.2 and 0.5. m^/mf is fixed to 1.0. 



As we are interested in the real part of the pole of mally introduce a cutoff A which will be taken to infinity 

we hereafter omit ir\ in the denominator. Each term in later, 
the integrand shows an ultraviolet divergence, and we for- 



h(P ,\P\) = 



(27T) 



lim 



dk 



IPI 



k In 



2u 



|P|fc 

mi 



Pn 



'2rn i, 



A*fc 



2// 



jp|fe 



Po 



P 2 

2 m b 



M6 



2k 2 



(A.2) 
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FIG. 15: Fermion dispersion curves. Left panel corresponds Nb/Nf = 0.1 and right one corresponds Nt/Nf = 1.0. (fcpa) 1 
are fixed to 0.5. And mb/mf — 1.0. 



The divergent term at A — > oo coming from the first term in the integrand is cancelled out by the second term, and 
we obtain the finite quantit 



I(Po-\P\) = 



2IPI 



m b 



A In 



HP| 



.4 



u\P\ 



A 



—^k F - — -AF(A) 



m b 



(A.3) 



where 



F(A) 



2k F 



(k F + ^-VA)(kr->^-VA) 



• arctan 



k F + - 



arctan 



k F — 



(A>0) 
(A = 0) 

(A<0) 



(A.4) 



with 



.4 



HP| 

m b 



2v ^P 



2m b 



(A.5) 



The final expression for T(P) is given in terms of I(P) 
by 



r(P) 



2na 



2na . 

1 I(Po,\? 

v 



(A.6) 
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